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Abstract 

A thermodynamic analysis of the harmonic oscihator is presented. Motivation for the study is 
provided by the blackbody radiation spectrum; when blackbody radiation is regarded as a system of 
noninteracting harmonic oscihator modes, the thermodynamics follows from that of the harmonic 
oscillators. Using the behavior of a harmonic oscillator thermodynamic system under an adiabatic 
change of oscillator frequency oo, we show that the thermodynamic functions can all be derived from 
a single function of uj/T, analogous to Wien's displacement theorem. The high- and low- frequency 
energy limits allow asymptotic energy forms involving T alone or uj alone, corresponding to energy 
equipartition and zero-point energy. It is suggested that the "smoothest possible" function which 
behaves as a — bx at small values of x and vanishes at large x is ae~^^/'^ because it is a monotonic 
analytic function every derivative of which is a multiple of the function itself. In this sense, it is 
noted that the Planck spectrum with zero-point radiation corresponds to the function satisfying 
the Wien displacement result which provides the smoothest possible interpolation between energy 
equipartition at low frequency and zero-point energy at high frequency. 
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I. INTRODUCTION 



Blackbody radiation holds an unusual place in the physics curriculum. Blackbody radi- 
ation provides a simple system which can be used as a model for thermodynamic analysis 
showing the derivation of the Stefan-Boltzmann law connecting total radiation energy and 
temperature [1]. Also, thermal radiation can be regarded as a collection of harmonic radi- 
ation modes[2] which can be treated by Doppler shift from a moving piston to obtain the 
Wien displacement law relating to the blackbody spectrum[3]. However, the full blackbody 
spectrum is regarded as inaccessible from thermodynamics; rather the blackbody spectrum 
reappears in statistical mechanics where it is used as an illustration of the breakdown of 
classical statistical mechanics and the need for energy quanta[4]. 

In this article we depart from the traditional points of view. We treat blackbody radia- 
tion by analyzing the thermodynamics of the harmonic oscillator. Although the harmonic 
oscillator is often discussed in statistical mechanics, it is not usually analyzed in thermody- 
namics[5]. We begin by noting the expression for the work done by a harmonic oscillator on 
slow change of its frequency. Then we use the laws of thermodynamics to give a simple but 
unfamiliar derivation of Wien's displacement law [6]. This law is thus seen to hold not only 
for electromagnetic radiation but also for any harmonic oscillator which is a thermodynamic 
system at temperature T. Next we discuss all the thermodynamic functions associated with 
a harmonic oscillator thermodynamic system. We note that the high- and low-temperature 
limits allowed by the Wien law correspond to zero-point energy and energy equipartition. 
This simple observation is again not traditional in the physics hterature. Finally, we suggest 
a natural interpolation between the high- and low- temperature forms allowed by the Wien 
displacement law. The simple interpolation leads to the Planck spectrum with zero-point 
energy. 

The unconventional treatment presented here provides a compact approach to the black- 
body problem through an elementary thermodynamic system which is simple enough for the 

textbooks. The analysis emphasizes the intimate connections between blackbody radiation 
and zero-point radiation, and also fits neatly with the modern realization that constant ac- 
celeration of a harmonic oscillator through zero-point radiation makes the oscillator energy 
equal to the Planck spectrum with zero-point radiation [7]. 

II. THERMODYNAMICS OF A HARMONIC OSCILLATOR 

A harmonic oscillator corresponds to a nonrelativistic point mass m moving in one di- 
mension in a potential V = (l/2)mcj^a;^. It can be described by the Lagrangian L{x,x) = 
(l/2)mi;^ — (l/2)ma;^a;^, and corresponds to harmonic oscillation x = {2U /mui)^^'^ cos{ujt+(j)) 
at energy Li and angular frequency uj. Under a slow change in the spring constant of the 
system corresponding to a change in the natural angular frequency u, the energy divided 
by the frequency lA/uj is a constant [8] and so the change in system energy is dlA = {U/uj)duj 
and the work dW done by the system on the external agency changing the frequency is 

dW = -{U/uj)duj. (1) 

Thus we regard X — —lA/uj as the generalized force associated with a change in frequency 

U). 

If a harmonic oscillator is weakly coupled to a heat bath at temperature T, then it will 
exchange energy with the heat bath. In equilibrium with the bath, the harmonic oscillator at 
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frequency 00 will have an average energy W (a;. T). Thus the oscillator forms a thermodynamic 
system suitable for treatment by the usual methods of thermodynamics involving (average) 
energy parameter a;, and associated work dW = —{U/u})dui. 



III. DERIVATION OF FUNCTIONAL DEPENDENCE CORRESPONDING TO 
WIEN'S DISPLACEMENT THEOREM 

When the oscillator is at thermal equilibrium in a bath at temperature T, the oscillator 
will have an (average) energy V{{uj.T) and an entropy S{uj,T) depending upon the temper- 
ature T and natural oscillation frequency u). Using the work expression in(l), the laws of 
thermodynamics for the oscillator give 

dQ = TdS{T, 00) = dU{T, 00) - {U/u)duj (2) 

This can be rewritten as 

Since the variables uj and T are independent, this requires 

Differentiating the first equation of (4) with respect to u and the second with respect to T 
gives 

^ d's _ d'u (ds\ d's _ d'u 1 (du\ 

dudT dudT' \dujj^ dTduj dTdu uj \dT J ^ ^' 

Subtracting the first equation of (5) from the second as as to eliminate the second-derivative 
terms, we have 

fdS\ 1 /dU\ 



Next using (6) together with the first equation of (4), we find 

'dS\ T fdS' 



du u \ dT 



(7) 



which has the general solution 

S{uj,T)=g{uj/T) (8) 

where g is an arbitrary function of the single variable uu/T. On the other hand if we use the 
second equation of (4), then we find from Eq.(6) 

which has the general solution 

U{u,T)^uf{u/T) (10) 

where / is an arbitrary function of the single variable lo/T. This corresponds to the infor- 
mation in the Wien displacement theorem. [3] Although the Wien theorem is often derived 
from the doppler shift of electromagnetic waves reflected from a moving piston in a cavity, 
our analysis shows that it holds in thermal equilibrium for any system described by a simple 
harmonic oscillator Lagrangian. 
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IV. THERMODYNAMIC FUNCTIONS FOR A HARMONIC OSCILLATOR 



The Wien displacement results S{u!,T) = g{uj/T) and U{u!,T) = ujf^uj/T) give con- 
straints upon all the thermodynamic functions for a harmonic oscillator. The thermody- 
namic potential[9] (j){uj/T), from which all the other thermodynamics functions may be 
derived, must be a function of the combination ou/T. The average oscillator energy U in 
thermal equilibrium follows as 

^(^^T)=T'(^) =-uj^\uj/T) (11) 

The entropy S of the oscillator is again a function olto/T, 

S{oj/T) = (P{lu/T) + U{uj, T)/T = (P{lu/T) - {uj/T)(P'{lu/T) (12) 
The Helmholtz free energy is directly related to the thermodynamic potential (f){u;/T) 

J^{uj,T) = -T(t){uj/T) (13) 
The generalized force X associated with a change in uj is 

n-m-T{?±)^ = ^'(.IT) (14) 
and the specific heat C at constant ui is given by 

C(a./r) = (§)^ = (15) 

Thus the equilibrium thermodynamics of a classical harmonic oscillator system is determined 
by one function, the unknown function (j){uj/T). 

V. WIEN DISPLACEMENT RESULT AND ZERO-POINT ENERGY 

There are two natural extremes for the oscillator energy given by the Wien displacement 
result in (11); one extreme makes the energy U{ui,T) independent of temperature T, and 
the other makes the energy lA{uj,T) independent of the natural frequency uo. 

We deal first with the temperature-independent energy. When the potential function 
(j)'{uj/T) = —const and so 

(j){u/T) = -const X (u/T) (16) 
then the oscillator energy in (11) takes the form 

U{uj, T) = Uzp{<jo) = const x u (17) 

This corresponds to temperature-independent zero-point energy. 

Wc note that substitution of the zero-point energy (17) into the first law of thermody- 
namics in the form 

dQ — dlA — {U/uj)dijj = d{const x cu) — {const x u!/u!)du! — (18) 



4 



tells us that no heat dQ enters the system on changing the natural frequency of the oscillator 
ui. Thus changes in zero-point energy occur without any change in the thermodynamic 
entropy S{u}/T) of the system. Indeed, we see that if is constant, then 02p(a;/T) = 
const X uj/T must be linear in its argument, and the entropy S in Eq.(12) vanishes for any 
function which is linear in its argument. 

Zero-point energy is random energy which is present even at zero temperature. Ther- 
modynamics allows the possibility of zero-point energy and experimental evidence (such as 
that for van der Waals forces) requires its existence. [10] It is natural to choose the unknown 
constant for the zero-point energy in (17) so as to fit the experimentally measured van der 
Waals forces. This corresponds to an oscillator energy 

U,p{uj) = {l/2)noj (19) 

where ^ is a constant which takes the value familiar for Planck's constant. 

VI. WIEN DISPLACEMENT RESULT AND ENERGY EQUIPARTITION 

The other extreme for the Wien displacement result (11) is the case where the oscillator 
energy depends upon the temperature but has no dependence upon the natural oscillator 
frequency u. Thus when (f)'{u;/T) — — const/ {u /T) in equation (11), then the oscillator 
energy is 

U{oj, T) = Urj{T) = const x T. (20) 

This is the familiar energy cquipartition law (proposed by Raylcigh and Jeans for low- 
frequency radiation modes) where the constant is chosen as Boltzmann's constant ks, 

Urj{T) = ksT. (21) 

In this case, an isothermal change of the natural oscillator frequency u produces no change 
in the oscillator internal energy. Rather, from (2), the isothermal work done on changing 
the natural frequency uj is provided by heat added which keeps the internal oscillator energy 
constant, 

dQ = TdS{uj/T) = dUnjiT) - {U/u)du;, constant T (22) 

Then 

dS^j = — {k:B/uj)doj, constant T, (23) 
and since we know the functional form S{uj/T), we have the familiar result 

Srj{u/T) = -ke \vl{u/T) + canst (24) 

Indeed if (t)'Rj{uj/T) = -A;s/(tu/T), then (t)Rj{uj/T) = -kBln{u/T) and the entropy in (12) 
takes the form (24). 

VII. USE OF NATURAL UNITS IN THE ANALYSIS 

In this paper we are not interested in the numerical evaluation of thermodynamic quan- 
tities but rather in the fundamental thermodynamic behavior. On this account we will 
measure all quantities in terms of energy and take the entropy as a pure number. Thus we 
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will take h—1 and measure frequencies in energy units. Also, we will take = 1 and mea- 
sure temperature in energy units. [11] Thus the limiting form corresponding to zero-point 

energy has 

<P,p{ij/T) = -{l/2){uj/T), U^piu) = {1/2)00 (25) 
while the limiting form corresponding to energy equipartition becomes 

Mlj/T) = - ln{Lj/T), Urj{T) - T (26) 



VIII. ASYMPTOTIC LIMITS FOR THERMAL OSCILLATOR ENERGY 

In general, the behavior of an oscillator system will depend upon both frequency uj and 
temperature T as in (11). In the limit as T ^ 0, we expect to recover the zero-point energy 
of the oscillator 

limT^oUiu, T) = limT-.Q[-uj^'{uj/T)] = U^p{u) = {l/2)u (27) 

and the associated thermodynamic potential 

0(cu/T) ^ 4>.p{^/T) = -(l/2)(cu/r) (28) 

In the limit a; — > 0, we expect to obtain the equipartition energy 

lim^^QU{uj, T) = lim^^o[-uj(t)'{uj/T)] = Urj{T) = T (29) 

and the associated thermodynamic potential 

<t>{uj/T) ^ (t>Rj{uj/T) = - ln(a;/r) (30) 

It is sometimes useful to make a distinction between the THERMAL energy Ut{uj, T) of 
an oscillator and the oscillator's TOTAL energy U{u,T). The thermal energy is just the 
(average) energy above the (average) zero-point energy 

Ut{uj, T) = U{u, T) - U,p{u) = -u(P'{u/T) - {l/2)uj = -u[ct>'{u/T) - 0l,(cc;/T)] (31) 

Although the total oscillator energy U{u!,T) is related to forces, it is only the thermal 
oscillator energy Ut{u!,T) which is related to changes in thermodynamic entropy since (as 
seen above) (f)zp{uj/T) does not give any thermodynamic entropy. 

IX. PLANCK SPECTRUM AS THE SMOOTHEST INTERPOLATION BE- 
TWEEN THE EQUIPARTITION AND ZERO-POINT LIMITS 

Although the thermodynamic forms given in Eqs.(ll)-(15) represent the information ob- 
tained from thermodynamic analysis, it is tempting to try to guess the full thermodynamic 
behavior which is chosen by nature. Now experiments indicate that nature indeed chooses 
the asymptotic forms given in Eqs.(27)-(30) corresponding to non-zero values of h and ks. 
Thus we ask whether there is any natural choice for behavior which connects these asymp- 
totic forms. One expects thermodynamics to involve smooth functions, and hence one might 
seek the "smoothest possible" interpolation between the extremes. 
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Is the idea of a "smoothest possible" interpolation well-defined? In general, this seems 
ambiguous. However, there are cases where the "smoothest possible" function appears 
absolutely clear. Thus if for small values of x an analytic function has the form f{x) = 
a — bx + 0{x'^) with a > 0, 6 > 0, and for large values of x the function f{x) 0, then the 
smoothest possible interpolation between the limits is f{x) = a exp{—bx / a) . This function 
meets the asymptotic limits and is a monotonic function whose nth derivative is just {—b/a)"' 
times the value of the function at that point. Thus every point x involves the same local 
functional behavior for f{x). Every point is equivalent. It is only modification of this 
functional form which introduces preferred values for x. 

Let us reconsider the asymptotic forms given in Eqs.(27)-(30). The energy forms are not 
useful since these do not give us asymptotic functions oicu/T. The limiting thermodynamic 
potentials 

Mij/T) = - ln(a;/r), cf>,p{u;/T) = -{l/2)uj/T (32) 

involve logarithmic behavior at the low- frequency limit. Since the logarithmic function is 
more complicated analytically than the exponential function, it is convenient to take the 
exponential of the negative of these functions and to consider 

exp[-0Rj(cc;/T)] = cu/T, exp[-0^p(o;/T)] = exp[c^/(2T)] (33) 

The exponentiation will not change the "smoothest possible" criterion required of the inter- 
polation. Thus according to (33) we are searching for the "smoothest possible" interpolation 
exp[— ^(cu/T)] between linear behavior uj/T at small argument and exponential behavior 
exp[a;/(2T')] at large argument 

exp[— 0(a;/T)] ^ a;/T for a;/T ^ 0, and exp[— 0(a;/T)] — exp[a;/(2T)] for tu/T — oo 

(34) 

We notice that the difference of the high-frequency limit and the desired interpolation has 
the asymptotic forms 

exp[a;/(2T)] - exp[-0(u;/T)] ^ 1 + a;/(2T) - u/T = 1 - uj/{2T) for u/T 



cxp[cj/(2T)] — exp[— 0(a;/T)] — >■ for u/T^oo 



(35) 



But this corresponds to exactly the case mentioned in the previous paragraph. Thus the 
"smoothest possible" interpolation is 



exp[uj/{2T)] - exp[-0(u;/T)] = exp[-u;/(2T)] 



(36) 



This imphes 



exp[-0(u;/r)] = exp[u;/{2T)] - exp[-u; / {2T)] = 2 sinh[u;/(2r)] 



(37) 



It is easy to see that the right-hand side of (37) has exactly the asymptotic forms de- 
manded in (33). Then taking the logarithm of (37) and using the thermodynamic relations 
given in (11) and (12), this smooth interpolation leads to the thermodynamic functions 



^Pzp 



In 



2 sinh 



2T 



(38) 
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) 



exp(cj/T) - 1 2 




u 1 



) 



(39) 



(40) 



Wc have labeled these thermodynamic functions with the subscript "Pzp" or "P" because 
they correspond exactly to the familiar Planck average oscillator energy including zero-point 
energy. As noted above, the entropy depends upon the Planck thermal spectrum Uri^^, T) 
but does not reflect the zero-point energy included in ZY(a;,T). 

X. COMMENTS ON THE INTERPOLATION FOR PLANCK'S SPECTRUM 

The analysis above shows that demanding the smoothest interpolation between the 
equipartition and zero-point limits suggests the Planck spectrum. This seems a surprisingly 
simple extrapolation from the thermodynamic analysis. One may wonder why such an ex- 
trapolation was not made a century ago. Indeed, many physicists are aware that Planck 
did arrive at the blackbody spectrum as an interpolation. However, Planck's interpolation 
did not come from the limits on the Wien displacement theorem. Rather, his interpolation 
involved a modification which combined earlier guesses at the high- and low-frequency parts 
of the thermal energy spectrum 1/(t{(jJ,T). [12] 

The simple interpolation made here was not made a century ago because physicists did 
not think in terms of a temperature-independent zero-point radiation. Indeed, even today, 
textbook discussions of blackbody radiation do not usually make any reference to zero- 
point radiation. [13] It is only much more recently, and in particular in connection with 
experimental measurements of Casimir forces, that physicists have taken seriously the zero- 
point energy suggested by the high-frequency limit of the Wien displacement law. 

Indeed the smooth extrapolation analysis for a harmonic oscillator given here fits neatly 
with the thermal effects of acceleration suggested by Davies and Unruh[7]. If one considers a 
charged harmonic oscillator undergoing constant acceleration through zero-point radiation, 
then the oscillator comes into equilibrium with the random radiation. When the oscillator 
has zero acceleration, its average energy U{u!,a = 0) depends only upon its frequency and is 
just zero-point energy U{uj,a = 0) = (1/2) hw. When the oscillator has a large acceleration 
a, then the average oscillator energy depends only upon its acceleration a and is independent 
of its natural frequency uj. Thus zero-point energy and energy equipartition form the natural 
extremes. And the function which connects these extremes is found to be exactly the Planck 
spectrum (39) with a temperature T = ha/2TTckB-[7] 

XI. CONCLUDING SUMMARY 

The harmonic oscillator provides a simple thermodynamics system involving energy U, 
temperature T, and oscillator frequency u. If one carries out a quasi- static change in 
the oscillator frequency lu, then thermodynamic analysis leads to results which correspond 
to the Wien displacement law. All of the thermodynamic behavior for the oscillator can 
be derived from a single thermodynamic potential function (^{cu/T) depending upon the 
single variable lo/T. If we consider the limits which make the oscillator energy U{u,T) 
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independent of one of its variables, then we find the extremes corresponding to zero-point 
energy and energy equipartition. Finally, as a guess for nature's choice of thermodynamic 
behavior, we ask for the function satisfying the Wien displacement law which provides the 
smoothest interpolation between the extremes of zero-point energy and energy equipartition. 
This leads to the Planck spectrum with zero-point radiation. 
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